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ABSTRACT 

We determine the metallicity distribution function (MDF) of the Galactic halo by means of a sample of 1638 metal-poor stars se- 
, lected from the Hamburg/ESO objective-prism survey (HES). The sample was corrected for minor biases introduced by the strategy 

for spectroscopic follow-up observations of the metal-poor candidates, namely "best and brightest stars first". Comparison of the 
metallicities [Fe/H] of the stars determined from moderate-resolution (i.e., R ~ 2000) follow-up spectra with results derived from 
abundance analyses based on high-resolution spectra (i.e., R > 20, 000) shows that the [Fe/H] estimates used for the determination of 
the halo MDF are accurate to within 0.3 dex, once highly C-rich stars are eliminated. We determined the selection function of the HES, 
which must be taken into account for a proper comparison between the HES MDF with MDFs of other stellar populations or those 
predicted by models of Galactic chemical evolution. The latter show a reasonable agreement with the overall shape of the HES MDF 
for [Fe/H] > -3.6, but only a model of Salvadori et al. (2007) with a critical metallicity for low-mass star formation of Z„ = 10 -3 ' 4 Z Q 

reproduces the sharp drop at [Fe/H] 3.6 present in the HES MDF. Although currently about ten stars at [Fe/H] < -3.6 are known, 

the evidence for the existence of a tail of the halo MDF extending to [Fe/H] ~ -5.5 is weak from the sample considered in this 
paper, because it only includes two stars [Fe/H] < -3.6. Therefore, a comparison with theoretical models has to await larger statis- 
tically complete and unbiased samples. A comparison of the MDF of Galactic globular clusters and of dSph satellites to the Galaxy 
shows qualitative agreement with the halo MDF, derived from the HES, once the selection function of the latter is included. However, 
statistical tests show that the differences between these are still highly significant. 
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1. Introduction 

One of the key observables for constraining models of the for- 
mation and chemical evolution of the Galaxy is the Metallicity 
Distribution Function (MDF) of the constituent stars of its vari- 
ous components (bulge, disk, halo). The MDF provides critical 
information on the enrichment history of those components with 
heavy elements. In the case of the halo, early enrichment may 
have been provided by the very first generations of massive stars, 
formed from material of primordial composition shortly after the 
Big Bang (i.e., Population III stars). 

Models of Galactic chemical evolution need to be compared 
to an accurate (and precise) observed halo MDF to test their pre- 
dictions, to constrain their various parameters (such as the ef- 
fective yield, the star-formation rate and the IMF), and in or- 
der to obtain information on the properties of Population III 
stars that are responsible for the earliest enrichment. This is par- 
ticularly important for the lowest metallicity tail of the MDF, 
which provide s invaluable in formation on the earliest enrich- 
ment phases (Prantzos 2003]); for instance, it has been sug- 
gested that a minimum level of enrichment is required to form 
low-mass stars. T his critical me t allicit y ranges between 10~ 4 Z F 



dOmukail l2000t iBr omm et aT] l200lt IBromm & Loebl [2 003. 
Umeda & Nomotol I20031 ISantoro & Shulll l2006t IFrebel etalj 



2007b and 1 6 Z , the latter being applicable when dust grains 



are present dSchneider et aT] 120021 12003L )2006t lOmukai et alj 
l200llTsuribe & Omukaill2006tlClark et al.ll2008l) . 

The precision of a derived halo MDF increases directly with 
the total number of observed metal-poor halo stars. Selection of 
such stars without the introduction of a kinematic bias (e.g., from 
among high proper motion stars) makes them of particular utility 
for examination of the relationships between the chemistry and 
kinematics of the halo. Early determinations of the halo MDF 
were based on small samples of globular clusters (lHartwickl 
119761: N = 60). o r a mixture of halo subdwarfs and globular clus- 
ters ([Bond 1981; N - 90 and N = 31, respectively). Problems 
with these samples arise not only from their small sizes, but 
also their inaccurate metallicities. Later studies employed signif- 
icantly larger samples wit h spectroscopically-de termined stellar 
abundances. For example. iRvan & Norris] ( | 1 99 lb used a s ample 
of 372 kinematically-se lected halo stars. IRvan & Norrisl d!991l) 
and lCarnev et al.l (Q~996) showed that the MDF peaks at a metal- 
licity of [Fe/H] = -1.6 with wings from [Fe/H] = -3.0 to solar 
abundances. 

The HK survey dBeers et al.lll985lfT99llBeersll 19991) . origi- 
nated by Preston and Shectman, and greatly extended by Beers to 
include several hundred additional objective-prism plates, was, 
until the advent of the Hamburg/ESO Survey (HES; see below), 
the primary source of metal-poor candidates suitable for con- 
sideration of the halo MDF. With the assistance of numerous 
colleagues, medium-resolution spectroscopy of over 10,000 HK- 
survey stars was obtained, using 1 .5^1 m class telescopes, over 
the past two decades. This led to the identification of thousands 
of stars with [Fe/H] < -2.0, as well as significant numbers of 
stars with [Fe/H] < -3.0. 

Another wide-angle spectroscopic survey is the HES. It was 
originally conceived as a s urvey for bright quasars dReimersI 
1990; IWisotzki et al.|[T996l 12000); however, its data quality is 
sufficient to not only efficiently select quasars with redshifts 
of up to z = 3.2, but als o various types of ste llar objects, in- 
cluding metal-poor stars (Christlieb et al. 2008). So far, several 
hundred new stars at [Fe/H] < -3.0 have been identified, in- 
cluding three stars that were confirmed by high-resolution spec- 
troscopy to have [Fe/H] < -4.0: HE 1327-2326 ([Fe/H] = 
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Fig. 1. Upper panel: Isochrones for an age of 12Gy r and metal- 
licities of [Fe/H] = -1,-2, and -3 dKim et al.l2002l) , and chosen 
colour cuts (see text for details); middle panel: V magnitude dis- 
tribution of the HES sample from which we construct the halo 
MDF; lower panel: (B - V)q distribution. 



5.4; iFrebel etal J 120051: lAoki et all 120061: iFrebel et alJl2006ah : 
HE 0107-5240 (fF e/H] = -5.3: IChristlieb et al J I2002T T2004: 
Bessell etall 120 04); and HE 0557-4840 ([Fe/H] = -4.8; 



Norri s et al. 1 120071 t is perhaps of interest that the HK survey 
has not (to date) yielded any stars with [Fe/H] < -4.0 confirmed 
by high-resolution spectroscopy; this may be related to the fact 
that the HK survey reaches apparent magnitudes that are brighter 
than the HES, and as a result is dominated more than the HES 

by inner-halo stars. 

The Sloan D igital Sky Survey (SDSS; iGunn et al.l I1998L 
lYork et all 120001) . and in particular the Sloan Extension for 
Galactic Understanding and Exploration (SEGUE) , has provided 
even larger samples of halo stars, as discussed by ICarollo et al.l 
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(2007j) and llvezic et al.l d2008l) . The former emphasize the divi- 
sion of the halo into two structural components, an inner region 
with R < 10-15 kpc, and an outer region beyond that radius. 
These two components differ in stellar metallicities, stellar or- 
bits, and spatial density profiles. As we discuss in Sect. |2]below, 
the HES sample is dominated by inner-halo stars. We note that 
we hereafter refer to the inner halo as "the halo", unless indi- 
cated otherwise. 

In spite of the very large sample of ~ 20, 000 stars used by 
Carollo et al., their coverage of the regime of very low metallic- 
ity is limited. According to their supplemental Fig. 4, they find 
only 3 stars with [Fe/H] < -3.0 in their "local sample" of 10,123 
stars. The main reason for this is that the stars of their sample 
were not selected to be metal-poor, but for the purpose of spec- 
trophotometric and telluric calibration of the SDSS spectra. 

Recent high-resolution spectroscopic follow-up of stars from 
the Carollo et al. sample (W. Aoki, priv. comm.) has indicated 
that the curr ent version of the S E GUE Stellar Parameter Pip eline 
(SSPP; see iLee et al.l l2008allbl: lAllende PrietoetaTI l2008h is 
somewhat conservative in the assignment of stellar metallicity 
estimates, in the sense that stars assigned [Fe/H] < -2.7 by the 
SSPP are in reality more metal-deficient, on average, by on the 
order of 0.3 dex. A recent examination of the numbers of stars 
from the SDSS/SEGUE survey, taking into account this offset, 
suggests that up to several hundred stars with [Fe/H] < -3.0 are 
in fact present in the current SDSS sample of stars (including 
other categories of targ ets than just the calibration stars). 

llvezic et all (12008) focus on the comparison between the in- 
ner halo and the disk. Since they rely on abundances determined 
from photometry, they cannot reliably determine metallicities of 
stars at [Fe/H] < -2. Nevertheless, the metallicity map of some 
2.5 million stars with photometric metallicies shown in Fig. 8 
of Ivezic et al. indicates that there exist very large numbers of 
stars in SDSS consistent with [Fe/H] < -2.0. Follow-up spec- 
troscopy is, at present, only available for a subset of them. Beers 
et al. (in preparation) discuss the MDF of the lowest metal- 
licity stars found in SDSS/SEGUE. The total number of stars 
with [Fe/H] < -2.0, based on medium-resolution SDSS spec- 
troscopy, is over 25,000 (i.e., five times the number discovered 
by the combination of the HK and HES). 

Th is paper continues our series on t he stellar content of the 
HES (IChristlieb et ail I2001U Paper I; IChristlieb et all 1200 lal 
Paper II; IChristliebe t al. 2005, Paper III; IChristlieb et al.ll2008l 
Paper IV). We are mainly concerned with the low-metallicity 
tail of the halo MDF, which is constructed from a sample of 
1638 metal-poor stars selected in the HES by quantitative cri- 
teria (Sect. [2J. The follow-up observations and determination of 
the metallicities are described in Sect. [3] In Sect. [4] we detail how 
the MDF was constructed. We discuss the shape of the halo MDF 
in Sect. [5] Comparisons of the observed MDF with MDFs pre- 
dicted by models of Galactic chemical evolution are presented in 
Sect. [6] and a comparison with the MDFs of the Galactic globu- 
lar cluster system and dwarf spheroidal galaxies is presented in 
Sect. |7] The results are discussed in Sect. [8] 

2. The metal-poor star sample 

One of the main advantages of the HES for determining the 
halo MDF is that the selection of candidate metal-poor stars 
was done with quantitative criteria. Hence, the selection is well- 
understood, and possible selection biases can be quantified and 
corrected for during the construction of the MDF. Futhermore, 
the selection is purely spectroscopic, so it does not introduce 
any kinematic biases. 



Table 1. Number of stars in each candidate class in the total sam- 
ple of candidates, number of observed candidates, and number 
of accepted candidates after removal of emission line objects, 
"peculiar" objects (e.g., objects with continuous spectra) and all 
stars with a G-band index GP > 6 A. In the last column, we list 
the scaling factors applied to the [Fe/H] histograms for each can- 
didate class during the construction of the MDF (see Section|4]i. 



Class 




Number of stars 


Factor 


All 


Observed 


Accepted 


mpca 


201 


123 


105 


1.63 


unid 


231 


208 


192 


1.11 


mpcb 


2006 


1008 


940 


1.99 


mpcc 


1275 


432 


401 


2.95 


Sum 


3713 


1771 


1638 





I 1 1 1 1 I 1 1 1 1 I 1 1 1 1 I 1 1 1 1 I 1 1 1 1 I 1 1 1 1 I 




10 20 30 40 50 60 

Distance [kpc] 



Fig. 2. Distance distribution of the HES sample. The sample is 
dominated by stars at distances of less than ~ 5 kpc from the 
Sun; a few cool giants are located at distances of up to ~ 50 kpc. 

The selection of candidates in the HES is described in 
Paper IV. For the sample used in this study, we employed only 
the KP/(B - V)n selection; i.e., a star is selected as a metal- 
poor candidate if its KP index of the Ca n K line, as measured 
in its digital HES objective-prism spectrum, is smaller than the 
KP index predicted for a star of [Fe/H] = -2.5 and the same 
(B - V)o colour (see Fig. 5 of Paper IV). This cutoff was cho- 
sen because it results in a good compromise between complete- 
ness at [Fe/H] < -3.0, the region in [Fe/H] we are mainly inter- 
ested in (because it corresponds to the earliest phases of Galactic 
chemical evolution), and achieving a selection that efficiently 
rejects stars at higher metallicity. In addition to the KP index, 
the B - V colours are measured in the HES spectra as well (see 
Paper IV for d etails), and then are c orrected for reddening us- 
ing the maps of Schlege Fet al.l (119981) . We restrict the sample to 
the colour range 0.5 < (B - V)o < 1.0, because the follow-up 
observations of stars bluer than (B - V)o = 0.5 have not yet 
reached a sufficient level of completeness, and for stars redder 
than (B - V)o = 1 .0, the accuracy of the determination of [Fe/H] 
from moderate-resolution follow-up spectra is limited due to the 
lack of calibration stars and the weakness of the H<5 line, which is 
used as a temperature indicator. The V magnitude and {B - V)o 
distribution of our sample together with isochrones for an age 
of 12Gyr and different metallicities is shown in Fig.[T] The V 
magnitudes as well as the (B - V)o colours are from the HES. 

The selection was applied to all spectra of unsaturated point 
sources extracted on 329 (out of 379) HES plates, covering a 
nominal area of ~ 7700 deg 2 of the southern high galactic lati- 
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Table 2. Follow-up observations of the candidate metal-poor stars. 



Telescope(s) 


Instrument(s) 


Observers 


Altars 


Magellan 1&2 


B&C 


Shectman, Mc William, Thompson 


553 


SS0 2.3m 


DBS 


Bessell, Norris, Edvardsson, 
Behnke, Christlieb, Frebel 


339 


Palomar 200" 


DS 


Cohen, Ramirez, Melendez 


323 


UK Schmidt 


6dF 


Haynes, Cass, Hartley, Russell, Watson 


283 


ES0 3.6m 


EFOSC2 


Fechner, Zickgraf, Barklem, 
Fuhrmeister, Christlieb 


140 


Total 






1638 
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Fig. 3. Spatial distribution of the HES sample. /?(gal. plane) is 
the distance from the Galactic center projected onto the Galactic 
plane. 



tude sky. The candidates were visually inspected and assigned 
to the classes mpca, unid, mpcb, and mpcc. As described in 
Paper IV, the classification is based on the appearance of the 
Can K line in the digital HES spectra. Candidates of class 
mpca are the best in terms of the success rate of finding stars 
at [Fe/H] < -2.5 (see Fig. |8j, since no Can K line could be 
seen in the HES spectrum, while the candidates of class mpcc 
are the worst, because a strong Ca K line could clearly be seen. 
However, the Ca K line is still strong in cool, moderately metal- 
poor (i.e., [Fe/H] 2.0) giants, therefore the line is expected 

to be detected in the HES spectra of such stars. For statistical 
studies such as the determination of the halo MDF it is therefore 
necessary to obtain follow-up spectroscopy also of the mpcc can- 
didates, because otherwise a color-related bias would be intro- 
duced. Furthermore, the assignment of classes to the candidates 



is subjective, and therefore it would be impossible to determine 
the selection function of the HES if only a subset of the candi- 
dates selected by quantitative criteria would be considered for 
the construction of the MDF. 

The result of the visual inspection are 3792 accepted can- 
didates, of which 79 are present on multiple plate quarters or 
plates; the number of unique candidates is 3713. The number 
of candidates in each of the aforementioned classes is listed in 
Tab. [TJ Only about half of the 3713 candidates are part of the 
sample presented in Tab. A. 1 of Paper IV. This is because slightly 
improved sky background and spectrum extraction algorithms 
were used in the final reduction of the HES, from which the 
sample of Paper IV was drawn. While minor changes of the re- 
duction algorithms can have a large effect on the measurement 
of the KP index in individual spectra, because the Ca n K line is 
covered by only four pixels of the HES spectra, we verified (HES 
plate by HES plate) that there are no systematic differences be- 
tween the KP indices measured in spectra reduced with the older 
extraction algorithms and the spectra to which the selection de- 
scribed in Paper IV was applied. Therefore, there should not be 
any statistical differences between the HES metal-poor sample 
presented in Paper IV and the sample used in this paper. We de- 
cided to construct the halo MDF from an older sample because 
the spectroscopic follow-up observations of that sample is more 
advanced, resulting in a considerably larger sample size. 

We determine distances to each of the sample stars using 
the [Fe/H] for each star and a set of isochrones similar to those 
shown in the upper panel of Fig.Q] Assuming that all the sample 
stars are at or above the main-sequence turnoff, we obtain the 
distance distribution shown in Fig. [2] and the spatial distribution 
shown in Fig. [3] The cooler giants in our sample reach distances 
from the Galactic plane well beyond |Z| = 15 kpc. However, the 
sample is clearly dominated by inner-halo stars. There is a hint 
that the outer-halo stars with |Z| > 15 kpc have a higher fraction 
of extremely metal-poor stars than do those of the inner halo 
with 5 < |Z| < 15 kpc, but given the wide range in metallicity we 
see throughout the halo, our sample is too small to determine the 
MDFs of the inner and outer halo separately with confidence. 



3. Determination of metallicities 

For 1771 of the 3713 unique candidates, moderate-resolution 
spectroscopy was obtained with various telescope/instrument 
combinations (see Tab.|2]l. The candidates were mostly observed 
in programs aiming at the identification of targets for high- 
resolution spectroscopy of the most metal-poor stars. Hence, the 
observing strategy adopted for the follow-up observations was to 
observe the brightest and best candidates (i.e., candidate classes 
mpca and unid) first. 
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Fig. 4. Pairs of KP and HP2 measurements for the same star in spectra obtained with different telescope/instrument combinations. 
Note that some of the estimates of <x displayed in the panels are influenced by a number of outliers; i.e., robust estimates would 
yield smaller values. 



In the follow-up spectra, we measured the KP index as well 
as th e HP2 index of U 6 and the GP index for the G-band of CH 
(see iBeers et al.lfl999l for the definition of these indices). When 
multiple spectra for a star were available, the S /W-weighted av- 
erage of the individual line index measurements was adopted. 
[Fe/H] was determined from the adopted HP2 and KP indices 



using the method of IBeers et al.l (1 19991) . Since the publication 
of that paper, the algorithm was improved mainly by including 
more calibration stars, which results in better coverage of the rel- 
evant stellar parameter space, and in particular in a better cover- 
age of the low- [Fe/H] region. 



6 



Schorck et al.: The halo metallicity distribution function 



£ -2 



b, -3 



-| i i i i | i i i i |- 
87 HERES stars 




-4 -3 -2 -1 

[Fe/H] high _ res 

Fig. 5. Comparison of determinations of [Fe/H] from moderate- 
resolution follow-up spectra with results based on high- 
resolution spe ctroscopy. Upper p anel: 87 stars observed with 
VLT/UVES dBarklem et all 12005). Lower panel: 23 stars ob- 
served with Keck/HI RES and two with Magellan/MIKE; anal - 
yses carried out b y ICohen et all d2004l) . ICohen et al] {2006), 
ICohen et a l. (2008), and Cohen (2008, unpublished). 



Since the stars of our sample were observed with many tele- 
scope/instrument combinations, it is important to verify that 
there are no systematic offsets between the measurements of 
the line indices in spectra taken at different telescopes. Such 
offsets could occur, for example, if the CCD response curves 
would strongly vary from instrument to instrument in the wave- 
length ranges in which the line and continuum bands of the 
indices are measured. For this reason, a number of candidate 
metal-poor stars were intentionally re-observed at different tele- 
scopes. Furthermore, in most of the observing campaigns, spec- 
tra of a few metal-poor standards (e.g., G 64-12, HD 140283, 
or CD -38° 245) as well as metal-poor radial-velocity standards 
were secured. In Fig. |H we show comparisons of the KP and 
HP2 indices measured in spectra taken with all relevant tele- 
scope/instrument combinations. In total, 315 pairs of measure- 
ments are available. No systematic offsets between the measure- 
ments can be seen. However, the scatter of the measurements in 
spectra obtained with the UK Schmidt and the fibre-fed multi- 
object spectrograph 6dF are about a factor two larger than those 
of the other telescope/instrument combinations. This can be at- 
tributed to the fact that sky subtraction is more difficult for the 
6dF spectra, since only a few fibers were dedicated to measure 
the sky background, and furthermore the sky brightness might 
have varied over the 6° diameter field of view of the instrument. 



The quality of the spectra (i.e., R ~ 2000 and a typical S IN 
of 20 per pixel in the continuum near the Ca K line) allowed the 
easy identification and rejection of emission-line and other "pe- 
culiar" objects (e.g., galaxies, or objects with continuous spectra, 
such as cool, heliu m-rich white dwarfs). It has been shown by 
ICohen et alj (120051) that CH lines present in the continuum bands 
of the KP and HP2 indices lead to a systematic underestima- 
tion of these indices, resulting in systematically too low [Fe/H] 
values. Hence, we also excluded from this study all stars with 
GP > 6 A. Since the fraction of carbon-enhanced stars among 
metal-poor stars increases as the metall icity decreases (see, e.g., 
ICohen et al] 120051: iLucatello et alj|2006l) . the rejection of stars 
with strong G-bands might lead to a bias against low-metallicity 
stars. However, since only 90 stars, or 5 % of the 1771 observed 
stars, were rejected due to this reason, the possible effect on 
our sample is only minor. We also note that the three currently- 
known ultra metal-poor stars (i.e., stars with [Fe/H] < -4.0; see 
Sect. |4] below), all of which have large overabundances of car- 
bon, are not rejected by this criterion, since their GP indices are 
smaller than 6 A. In total, 133 stars were rejected, leaving 1638. 

Homogeneous abundance analyses based on high-resolution 
spectra are available for 112 of the confirmed candidates in 
our sample. The spectra were taken with VLT/UVES (87 
stars), Keck/HIRES (23 stars) or Magell an/MIKE (2 sta r s). Th e 



abun dance analyses were performed b y Barklem et al 



ICohen et all d2004l) . ICohen et all d2006l) . ICohen et al. 



(2005), 
(2008) 



and Cohen (unpublished). Fig. compares the iron abundances 
determined in the course of these analyses ([Fe/H] high _ res ) to the 
moderate-resolution follow-up results ([Fe/H] med _ res ). No sig- 
nificant trends or offsets are present, and the 1-cr scatter around 
a regression line of the combined test sample is 0.3 dex. We 
hence conclude that the accuracy of [Fe/H] med _ res for our sam- 
ple is +0.3 dex. We note that the accuracy can be increased espe- 
cially for the cooler stars by using CCD photometry rather than 
B - V colors predicted from the HS index HP2 when deriving 
[Fe/H] med _ res . However, CCD photometry is not yet available 
for all stars of our sample. 

To increase the accuracy of the determination of the shape of 
the low-metallicity tail of the MDF, we replaced [Fe/H] med _ res 
with [Fe/H] high _ res , where available. [Fe/H] high _ res values are 
available for 27 of the 76 stars at [Fe/H] med _ res < -3.0, and five 
out of the six with [Fe/H] med _ res < -3.5. The [Fe/ H] higll _ rt , s val- 
ues we re taken from the references above and from Cavre l et al.l 
(120041) for HE 0305-5442, a re-discovery of CS 22968-014 
([Fe/H] high _ res = -3.56). The sixth star at [Fe/H] med _ res < -3.5 
in our sample for which a [Fe/H] high _ res estimate is available has 
[Fe/H] med _ res = -4.2. A VLT/UVES spectrum exists for this 
star, and a preliminary abundance analysis confirms that the star 
has a metallicity close to or slightly below [Fe/H] = -4.0. Due 
to the preliminary nature of this result, we do not show this data 
point in Fig. [5] 



4. Construction of the observed MDF 

In order to investigate potential selection biases given the 
adopted follow-up observation strategy, it is instructive to com- 
pare the MDFs derived from stars of the individual candidate 
classes and in different magnitude ranges. For the purpose of in- 
vestigating the possible presence of a bias caused by the fact that 
the brightest stars were observed first, we divided the full HES 
sample, as well as the subsamples of the four candidate classes, 
into a bright (B < 16.7) and a faint (B > 16.7) half, respectively. 
The results are shown in Fig. [7] 
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ity demonstrate that the candidate classification effectively sep- 
arated the "good" from "bad" candidates. Fig.[8]also shows that 
the number of false positives (i.e., stars at [Fe/H] > -2.5) is 
considerably higher among the mpcc candidates. However, this 
contamination does not affect our study, because we are mainly 
concerned with the low-metallicity tail of the MDF. 



o 1 ■ 1 1 ■ 1 1 1 ■ ■ 1 ■ ■ ■ ■ 1 ■ ■ ■ ■ 1 1 ■ ■ ■ 1 ■ ■ ■ ■ 1 ■ ■ ■ ■ 1 1 i 1 ■ 1 ■ ■ ■ ■ 1 ■ ■ ■ 1 1 

13 13.5 14 14.5 15 15.5 16 16.5 17 17.5 18 
B HES 

Fig. 6. Fraction of the stars in the HES sample for which 
moderate-resolution follow-up spectroscopy exists as function 
of B magnitude. 



In the upper right panel it can be seen that the faint can- 
didates are over-represented in the class unid. The reason is 
that the visual classification for fainter candidates, which have 
lower quality HES spectra, was more difficult. The bright- and 
faint-star [Fe/H] distributions of the other candidate classes, as 
well as the total sample, appear very similar to one another. This 
is quantitatively confirmed for the mpca and unid subsets by 
means of a Kolmogorov-Smirnov (KS) test of the null hypoth- 
esis Ho that the bright and faint subsets of the stars belonging 
to these candidate classes were drawn from the same parent dis- 
tribution, against the alternative hypothesis Hi that they were 
not drawn from the same parent distribution. We chose a signifi- 
cance level of a = 0.01; i.e., Ho is rejected if the probability p of 
occurence of the test statistic (i.e., in case of the KS-test, D, the 
maximum distance between the cumulative probability distribu- 
tions of the two samples), given Ho, is smaller than 0.01. The 
result of the KS-test for the bright and faint stars of the classes 
mpca and unid are p — 0.82 and p = 0.21, respectively; i.e., 
Ho can clearly not be rejected in these cases. However, for the 
other two classes, the probabilities are considerably lower, rang- 
ing from 0.0046 (mpcb) to 0.073 (mpcc). For the full sample 
(i.e., all candidate classes combined), the probability is 0.0028; 
that is, the distributions of the bright and faint subsamples differ 
significantly from each other. It would hence be desirable to con- 
struct the halo MDF from the bright and faint samples separately. 
However, in this case the sample sizes would be too small to 
draw any conclusions, in particular about the low-metallicity tail 
of the MDF. Therefore, we decided to construct the halo MDF 
from the full sample. We note that the relative fraction of ob- 
served stars does not exhibit any strong biases towards bright or 
faint stars (see Fig. [6j, thus the sample from which we construct 
the MDF should at least be representative for halo stars in the 
magnitude range 13 < B < 17.5. 

As can be seen in Fig. [8] the fraction of stars at [Fe/H] < -3.0 
is highest among the mpca candidates (i.e., 7 %), and signifi- 
cantly lower in the other classes (i.e., 3^4 %). That is, the MDF 
of mpca candidates is biased towards lower metallicities. KS- 
tests show that the null hypotheses Ho that subsamples of dif- 
ferent candidate classes were drawn from the same parent distri- 
bution can clearly not be rejected for neighbouring classes (e.g., 
p = 0.34 for unid versus mpcb), while Ho can be rejected at high 
significance when more distant classes are compared to each 
other (e.g., p — 1 .7 ■ 10~ 5 for mpca versus mpcc). These tests and 
the bias of the candidates of class mpca towards low metallic- 




Fig. 8. Metallicity distribution of the HES sample of 1638 stars, 
divided by candidate class. In the upper left corner of each panel, 
the number of stars with [Fe/H] < -3.0 and the percentage of 
such stars within each candidate class is indicated. 
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Fig. 7. Comparison of the MDFs of the bright (B < 16.7) and faint (B > 16.7) subsamples for each of the four candidate classes 
(upper four panels) as well as for the combined candidate sample (lower left panel). At the top left of each panel, the number of 
candidates belonging to the bright and faint sample, respectively, is listed; at the top right, the mean [Fe/H] of the samples is given. 

Table 3. The MDF of the Galactic halo field stars as constructed from the sample of 1638 HES with available spectroscopic 
follow-up observations, by means of scaling to the full candidate sample of 3439 stars (for details see text). Note that for a proper 
comparison with the MDFs predicted by theoretical models, or the MDFs of other stellar populations, the selection efficiency of the 
HES as a function of [Fe/H] and (B - V)q must be taken into account (see Tab.|4|. Note also that the HES sample is increasingly 
contaminated with thin- and thick-disk stars as [Fe/H] approaches 0.0. 
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In order to properly take into account the stars of our candi- 
date sample for which no spectroscopic follow-up observations 
exist, we constructed MDFs from the observed sample of stars 
in the following two ways. First, we computed separate MDFs 
for each of the candidate classes and scaled them such that the 
correct relative fraction of stars is produced when the four scaled 
MDFs are coadded; i.e., the scaling factors listed in the last col- 
umn of Tab. Q] were applied. Secondly, we assigned to each of the 
1942 stars in the full candidate sample lacking follow-up obser- 
vations the [Fe/H] of a randomly selected star of the same candi- 
date class for which a follow-up spectrum is available. We also 
randomly rejected stars with a too strong G-band and "peculiar" 
stars according to the probabilities determined from the sample 
for which follow-up observations exist. In this way, a sample of 
3439 stars with the correct relative fraction of the candidates of 
the four classes was created. 
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o.oi 



1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 

scaled (factor) 
: — scaled (random) 



1 I 1 1 1 1 I 1 1 1 1 I 1 1 
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L±J 



L±J 



-4.5 



-3.5 



-3 -2.5 -2 
[Fe/H] 



-1.5 



-0.5 



Fig. 9. Comparison of the MDFs constructed from the HES sam- 
ple by means of random scaling and co-addition of the class-wise 
MDFs (solid black line) and scaling by factors (grey dotted line). 



The MDFs produced by these two methods are expected to 
be very similar to each other, because in each of them, the class- 
wise MDFs are scaled and then added to produce the final MDF; 
only the scaling methods are slightly different. Indeed, as can be 
seen in Fig. |9j the results do not differ significantly from each 
other. A^- 2 -test of the null hypotheses Ho that the two samples 
are drawn from the same distribution yields a probability of p = 
1.0; i.e., Ho can very clearly not be rejected. We adopt the MDF 
constructed by means of scaling the class-wise MDFs by a factor 
and adding them up. For this MDF, the numbers of stars in each 
metallicity bin are listed in Tab. [3] 

Note that we have not corrected the MDF for the fact that 
as metallicity decreases, given that the giants become brighter 
and the dwarfs fainter (see Fig.[T]i, the relative volumes surveyed 
in our magnitude limited sample become larger and smaller, re- 
spectively. At (B - V)p = 0.6, for e xample, the data of the Yale- 
Yonsei isochrones dKim etal.ll2002h for an age of 12Gyr sug- 
gest that the ratios of volume surveyed at [Fe/H] = -3 to that at 
[Fe/H] = -2 are 3.0 and 0.67 for giants and dwarfs, respectively. 
Due to the very small survey volume for dwarfs, no cool main- 
sequence star has so far been identified in the HES, and the sam- 
ple considered in this paper is dominated by giants. Therefore, a 
survey volume correction would lead to a reduced relative num- 
ber of stars at the lowest metallicities. 
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Fig. 10. Comparison of the MDF constructed from the HES 
sample (solid l ine) w ith that of the HERES sample analysed by 
Barklem et al. (2005, grey dotted line). The latter sample is bi- 
ased against stars at [Fe/H] > -3.0, because most stars at higher 
metallicity were intentionally removed due to the science aims 
of that survey. 



5. The shape of the halo MDF 

A prominent feature in both of the scaled MDFs is a sharp drop 
at [Fe/H] ~ -3.6 (see Fig. [9); in our (scaled) sample, only 
two out of 3439 st ars hav e rFe/H ] < -3.6. Such a drop was 
also recognized by iNorrisI (Q~999), and it has been seen in the 
Hamburg/ESO R-process En hanced star Survey (HERES; see 
Fig. 2 of iBarklem etalj|2005l and our Fig. [10l>. It reflects the 
fact that only very few stars at [Fe/H] < -3.6 were found in 
projects aiming at the identification and detailed study of the 
lowest metallicity stars of the Galactic halo, despite the consid- 
erable effort expended to find them (see, e.g. JCohen e t al. 2008, 
and references therein). 
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Fig. 11. Comparison of the hal o MDF constructed from the HES 
sample (histogram) with that of R yan & Norris! d!991l) . scaled to 
match the HES MDF in the range -3.4 < [Fe/H] < -2.5 (dotted 
line). 



The shape of the low-metallic ity end of the halo M DF could 
not be determined precisely by iRvan & Norris! d!991l) due the 
limited size of their sample, which contains only four stars at 
[Fe/H] < -3.4, and none with [Fe/H] < -4.0. As can be seen 
in Fig. Q3] in the range -3.4 < [Fe/H] < -2.5 their halo MDF 
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agrees very well with the HES MDF. In Fig. QT| one can see 
a disagreement between the two MDFs in the bin centered on 
[Fe/H] = -3.5; i.e., the number of stars at this metallicity in the 
sample of Ryan & Norris is higher by about a factor of two as 
compared to the HES sample. Alternatively, the number of stars 
in the range -3.4 < [Fe/H] < -2.5 (i.e., the metallicity range 
which has been used to scale the two MDFs onto each other) 
are under-represented in the sample of Ryan & Norris, or over- 
represented in the HES sample. Even though the number of stars 
at [Fe/H] < -3.4 in both samples is small, the difference is sig- 
nificant. A KS-test of the null hypothesis Hq that the HES sample 
and the sample of Ryan & Norris have the same parent distribu- 
tion at [Fe/H] < -2.5 yields a probability of p = 0.0087; i.e., 

must (barely) be rejected if a = 0.01 is chosen. The proba- 
bility increases to 0.0091 if the two stars at [Fe/H] < -4.0 are 
excluded from the HES sample. The reason for the discrepancy 
is currently unclear, but one might speculate that it is related to 
the kinematic selection of the sample of Ryan & Norris and/or 
a larger fraction of stars belonging to the outer halo population 
being present in the HES sample. 

Another feature of the halo MDF is a lightly populated 
tail extending to [Fe/H] < -5.0. The evidence for this fea- 
ture from our (scaled) sample alone is weak, since it con- 
tains only two stars at [Fe/H] < -3.6, and none at [Fe/H] < 
-4.3. However, currently some ten stars with [Fe/H] < -3.6 
have published abundance an alyses based on high-reso lution 
spectroscopy (see Tab. 4 of iBeers & ChristliebTl2005l for a 
recent review), including three additional stars at [Fe/H] < 
-4.0: HE 1327-2326 ([Fe/H] = -5.4; iFrebel et al.1 120051: 
lAoki et al.l2006llFrebel et al.ll2006al) , HE 0107-524 ([Fe/H] = 
-5.3: IChristlieb et alJ 120021 12004 feessell et al I 120041) . and 
HE 0557-4840 ([Fe/H] = -4.8: lNorris e"t"alll2007l) . These three 
stars are not part of our sample due to a variety of reasons. 
HE 1327-2326 is part of the bright HES me tal-poor sample con - 
sisting of stars above a saturation threshold (Frebel_etaL|2006b), 
while only unsaturated point sources entered the sample of this 
work. HE 0107-5240 was selected in a previous version of the 
candidate selection which was slightly less restrictive than the 
one we use here; as a result, this star misses the selection cutoff 
of 3.9 A for its HES {B - V) colour of 0.6 mag by 0. 1 A (i.e., the 
KP index measured in the HES spectrum is 4.0 A). And finally, 
HE 0557-4840 is located on one of the 50 HES plates which are 
not considered here. In conclusion, for an accurate determina- 
tion of the shape of the MDF at [Fe/H] < -4.0 it is required to 
compile even larger statistically complete samples of metal-poor 
stars. 



6. Comparison between theoretical MDFs and the 
halo MDF 

In a comparison of the observed MDF with MDFs predicted by 
theoretical models, one has to take into account the modification 
of the shape of the MDF by the selection of metal-poor candi- 
dates employed in the HES. In particular, uncertainties cr^ and 
ctb-v of the measurements of the KP index and B - V in the HES 
spectra result in a scatter of stars with [Fe/H] > -2.5 into the 
sample, and stars with [Fe/H] < -2.5 out of the sample. Each 
theoretical MDF under investigation is therefore converted into 
an MDF as it would be observed in the HES, by applying the 
metal-poor star selection criteria used in the HES. 

The first step in the conversion of a theoretical MDF is the 
simulation of a sample of stars with a distribution in [Fe/H] ac- 
cording to that of the theoretical MDF under investigation. The 
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Fig. 12. Selection function for HES metal-poor candidates of 
(B-V)o = 0.5, 0.7, and 1 .0, as determined from a simulated sam- 
ple of stars following the MD F predicted by the Simpl e Model of 
Galac tic chemical evolution (Searle & Sar genl|[l972l ; lHartwickl 
[19761) . 



[Fe/H] values are then conv erted into pa i rs of K P and (B-V)o by 
inverting the calibrations of lBeers eta El (119991) . Then, a subsam- 
ple was selected such that it follows the distribution in (B-V)o of 
the HES sample (see Fig.Q]). Taking into account the distribution 
in (B - V)o is important because the shape of the selection func- 
tion is determined by <x KP , <x B _y, and the gradient of [Fe/H] in 
the KP versus (B - V\) parameter space (see Fig. 4 of Paper IV); 
it varies with (B - V)o, as can be seen in Fig.[T2l 

The reader will note that we have excluded stars with 
(B - V)u < 0.5 from our sample, which in principle will af- 
fect the relative proportions of dwarfs admitted to our sample 
as a function of [Fe/H]. In practice, however, this is not a se- 
rious effect if we restrict our abundance range to abundances 
[Fe/H] < -2.0. Consideration of the Yale-Yonsei isochrones for 
an age of 12 Gyr, and for the Salpeter mass function (x = 1 .35), 
shows that the percentages of dwarfs with (B - V)o < 0.5 rel- 
ative to all main-sequence stars with mass greater that 0.40 M Q 
are 4, 13, 19, 22, 24, and 24 for [Fe/H] = -1.0, -1,5, -2.0, 
-2.5, -3.0, and -3.5, respectively. That is to say, the proportion 
of excluded dwarfs is relatively constant for [Fe/H] < -2.0. 

The next step in the procedure of converting theoretical 
MDFs into an MDF as it would be observed in the HES was to 
add random Gaussian errors with standard deviations according 
to the known measurement uncertainties <x KP , ctb-v in the HES 
to KP and (B - V)o assigned to each star. Finally, the KP/(B - V)q 
selection criterion was applied to the simulated sample of stars. 
The [Fe/H] distribution of the selected stars is the MDF as it 
would be observed in the HES. 

For the convenience of the reader, we list in Tab. [4] and show 
in Fig. Q~2] the HES metal-poor star selection function as de- 
termined with a simulated sample of stars following the MDF 
predicted by the Simple Model of Gala ctic chemical evolution 
dSearle & Sargentlll972l; lHartwickl 1 19761) . That model assumes 
that a fiducial "closed box" of primordial gas is enriched by suc- 
cessive stellar generations. Further model assumptions are that 
(i) the gas is well-mixed at all times (i.e there is a unique age- 
metallicity relation for the stars formed from that gas) and (ii) the 
stellar initial mass function (IMF) does not change with time. 
Analytical solutions can only be obtained if it is assumed that 
the evolutionary timescales of the enriching stars are neglible 
(the so-called Instantaneous Recycling Approximation or IRA). 
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Table 4. Selection function for HES metal-poor candidates in 
the colour range 0.5 < (B - V)o < 1 .0, as determined from a 
simulated sample of stars following the MDF p redicted by the 
Simp l e Model of Ga lactic chemical evolution dSearle & Sa rgent 
U972tlHartwicklll976l) . 



[Fe/H] 




Selected fraction at (B - 


V)o 




0.5 


0.6 0.7 0.8 


0.9 


1.0 



-4.05 
-3.95 
-3.85 
-3.75 
-3.65 
-3.55 
-3.45 
-3.35 
-3.25 
-3.15 
-3.05 
-2.95 
-2.85 
-2.75 
-2.65 
-2.55 
-2.45 
-2.35 
-2.25 
-2.15 
-2.05 
-1.95 
-1.85 
-1.75 
-1.65 
-1.55 
-1.45 
-1.35 
-1.25 
-1.15 
-1.05 
-0.95 
-0.85 



1.000 
0.958 
1.000 
0.961 
0.982 
0.954 
0.920 
0.924 
0.901 
0.861 
0.816 
0.744 
0.668 
0.563 
0.455 
0.340 
0.232 
0.140 
0.075 
0.034 
0.013 
0.004 
0.001 
0.000 
0.000 
0.000 
0.000 
0.000 
0.000 
0.000 
0.000 
0.000 
0.000 



1.000 
1.000 
0.962 
1.000 
0.987 
0.991 
0.991 
0.979 
0.976 
0.950 
0.919 
0.869 
0.801 
0.700 
0.583 
0.457 
0.337 
0.234 
0.149 
0.088 
0.046 
0.021 
0.009 
0.003 
0.001 
0.000 
0.000 
0.000 
0.000 
0.000 
0.000 
0.000 
0.000 



1.000 
1.000 
1.000 
1.000 
1.000 
1.000 
0.997 
0.993 
0.989 
0.984 
0.958 
0.928 
0.879 
0.812 
0.715 
0.592 
0.462 
0.331 
0.222 
0.136 
0.079 
0.043 
0.022 
0.011 
0.006 
0.003 
0.001 
0.001 
0.001 
0.000 
0.000 
0.000 
0.000 



1.000 
1.000 
1.000 
1.000 
1.000 
1.000 
0.996 
0.996 
0.991 
0.974 
0.954 
0.908 
0.852 
0.768 
0.658 
0.537 
0.406 
0.297 
0.203 
0.130 
0.081 
0.048 
0.028 
0.017 
0.010 
0.006 
0.004 
0.003 
0.002 
0.001 
0.001 
0.000 
0.000 



1.000 
1.000 
1.000 
1.000 
1.000 
1.000 
1.000 
0.997 
0.985 
0.981 
0.953 
0.900 
0.839 
0.743 
0.617 
0.488 
0.364 
0.264 
0.187 
0.127 
0.085 
0.058 
0.041 
0.028 
0.020 
0.015 
0.010 
0.007 
0.004 
0.002 
0.001 
0.000 
0.000 



1.000 
1.000 
1.000 
1.000 
1.000 
1.000 
1.000 
1.000 
1.000 
0.997 
0.983 
0.949 
0.918 
0.822 
0.709 
0.573 
0.433 
0.318 
0.217 
0.152 
0.100 
0.070 
0.051 
0.036 
0.028 
0.021 
0.017 
0.014 
0.008 
0.008 
0.004 
0.000 
0.000 



Such solutions can be generically obtained in the case of a closed 
box, and in some particular cases of outflow (gas loss from the 
box) and infall (gas flows into the box). Since the IRA turns out 
to be a very good approximation for elements ejected by massive 
stars, those analytical solutions can provide a powerful tool for 
the study of Galactic systems. 

In the framework of the Simple Model, the shape of the MDF 
can be described in terms of a unique parameter, the "yield", 
which is the ratio of the mass of newly-created metals to the 
mass locked in long-lived stars and stellar remnants. This is a 
very useful parametrization, because it is independent of the star 
formation history of the system (the major unknown in Galactic 
evolution studies). In the closed box model the yield depends 
only on the IMF (referred to as the "true yield"), while in the 
case of gaseous flows (infall and outflow) it depends also on their 
magnitude; this "effective yield", y e ff, is always smaller than the 
true yield. It turns out that the MDF peaks at a metallicity equal 
to the effective yield; this simple result allows one to determine 
the effective yield and to constrain the underlying physics (IMF, 
outflow rate, etc.) 

In Fig. [13] we compare the MDF of a Simple Model with 
y eS = -1.7 with the MDF observed in the HES. The HES MDF 
shows an excess of stars in the range -3.5 < [Fe/H] < -3.0. 
Alternatively, if the MDF of the Simple Model would be scaled 
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Fig. 13. Comparison of the MDF of a Simple Model with 
y e ff = -1.7 (black line) with the MDF constructed from the HES 
sample (histogram). The grey line illustrates how the MDF of the 
Simple Model would be observed in the HES. For the compari- 
son between the theoretical model and the observations, only the 
grey line in the metallicity range below [Fe/H] ~ -2.0 should be 
taken into account, because at higher [Fe/H] the HES sample is 
contaminated with thin- and thick-disk stars. 



such that it matches the observed MDF in this range, a large 
deficit of the number of observed stars in the range -3.0 < 
[Fe/H] < -2.0 with respect to the Simple Model would result. It 
is also neither possible to reproduce with the Simple Model the 
sharp drop of the observed MDF at [Fe/H] = -3.6, nor the tail 
at [ Fe/H] < - 3.6. 

iPrantzosI (120031) developed a modification of the Simple 
Model, which includes early infa ll, and later outfl ow of gas; the 
IRA is also relaxed in his model. IPrantzosI (12007) suggested that 
since the halo of the Galaxy has been assembled by merging of a 
large number of fragments, the MDF of the Galactic halo can be 
seen as the sum of the MDFs of these fragments. In his model, 
the chemical evolution histories of each of the fragments are 
still described by the Simple Model, using the observed mass- 
metallicity relation of dwarf galaxies to derive individual effec- 
tive yields. The halo MDF is then produced by integrating over 
a mass function of the fragments dete rmined in numerical sim- 
ulations. The MDFs of the models of dPrantz os 2003, 120071) are 
shown in Fig. Q3] Both MDFs match the HES MDF well in the 
range -3.5 < [Fe/H] < -2.5 and at [Fe/H] ~ -4.0, but the sharp 
drop at [Fe/H] 3.6 is not predicted by them. 

The next set of models that we consider are those of 
ISalvadori et al.l (120071) . who reconstruct the merger tree of the 
Milky Way with a semi-analytic approach including a chemi- 
cal evolution code. A free parameter in this model is the critical 
metallicity for low-mass star formation, Z cr . As can be seen in 
Fig- [H] the model with Z cr = 1O~ 34 Z reproduces the drop of 

the observed MDF at [Fe/H] 3.6 rather well. However, the 

model predicts that no stars at [Fe/H] < -3.6 should exist, while 
there are two such stars in our sample, and for about ten addi- 
tional stars in this metallicity range abundance analyses based 
on high-resoluti on spectroscopy have bee n published (see Tab. 4 
of the review of lBeers & Christliebll2005l) . 

The Salvadori et al. model with Z cr = 1O~ 4 O Z matches the 

HES MDF at [Fe/H] 4.0, but disagrees in the range -4.0 < 

[Fe/H] < -3.6, where ~ 30-70 stars are predicted, but none are 
present in our sample. The model with Z cr = over-predicts the 
number of stars in this metallicity range even more strongly, and 
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Fig. 14. Comparison of the MDF constructed from the HES 
sample (histogram) with models of ISalvadori et al.l (|2007) with 
different critical metallicities Z cr (black and grey lines). The 
region between the dashed lines indicates the uncertainty of 
the models due to different hierarchical merger histories of the 
Galaxy. 



it greatly over-predic ts that number of stars at [Fe/H] < -4.0, as 
already discussed by ISalvadori et al.1 (120071) . 

Finally, we compare in Fig.[T6]the HES MDF wit h that pre- 
dicted by the stochastic chemical enrichment model of lKarlssonl 
d2006l) . While the model matches the HES MDF at [Fe/H] ~ 
-4.0, and predicts the tail at [Fe/H] < -4.0 that is known to 
exist from additional stars published in the recent literature, the 

drop of the observed MDF at [Fe/H] 3.6 is not present in the 

theoretical MDF 
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Fig. 15. Comparison of the MDF constructed from the HES 
sample (histogram) w ith theore t ical pr edictions (black and grey 
lines). Upper panel: [PrantzosJ (!2003): lower panel: iPra ntzos 
(120071) . 



1000 - 



100 



10 



1 - 



I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I 
Karlsson (2006) 



Karlsson (2006) as observed in HES 
HES sample 




: ' ' I ' ' ' ' I ' ' I I I ' ' ' I I ' ' ' ' I ' ' ' ' I ' ' ' ' I ' ' ' ' I ' ' ' ' I ' ' ' ' I ' ' ' ' I ' 
-5 -4.5 -4 -3.5 -3 -2.5 -2 -1.5 -1 -0.5 
[Fe/H] 

Fig. 16. Comparison of the halo MDF constructed from the HES 
sample with the M DF predicted by the stochastic chemical en- 
richment model o f lKarlssonld2006h . 



To quantify our comparisons of the HES MDF with those 
predicted by the theoretical models discussed above, we car- 
ried out KS-tests of the null hypotheses Hq that the HES sam- 
ple and the individual samples simulated according to the MDFs 
of the models were drawn from the same parent distribution, 
at a significance level of a — 0.01. The tests were restricted 
to [Fe/H] < -3.0, because we are mainly concerned with the 
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shape of the low-metallicity tail of the MDF. The result of the 
tests are that Ho can not be rejected (p - 0.063) only in the case 
of the Salvadori et al. model f or Z cr = 10~ 3 4 Z o . However, we 
n ote that the statisti cal test of iKuiper (1962), which according 
to lPress et al.l (1 19921) is more sensitive than the KS-test to differ- 
ences at the ends of the two distributions under comparison (i.e., 
at the lowest and highest metallicities), yields p = 6.5 ■ 10~ 4 ; 
i.e., if this test is employed, Hq would be rejected at high signifi- 
cance. All other models considered yielded p < 10 -3 , regardless 
of which of the two tests were applied. 

7. Comparison of the halo field star MDF with that 
of other stellar populations 

It is of great interest to compare the halo MDF with the MDF 
found for other stellar populations, in particular for the system 
of Galactic globular clusters (hereafter GCs) and for the stars 
in dwarf spheroidal (dSph) galaxies. Since the most metal-poor 

Galactic GC has [Fe/H] 2.5, we need to establish whether or 

not there is a real deficit of GCs at lower Fe-metallicities com- 
pared to the halo field. 

For a proper comparison of the HES MDF with that of other 
stellar populations, it is mandatory that the selection function of 
the HES, as listed in Tab. [4] be taken into account. The values 
in that table can be used to correct the observed MDF for the 
selection of metal -poor candidates employed in the HES. This 
is particularly important at [Fe/H] > -2.5, where the corrections 
are large, because typically less than half of the stars are actually 
picked up by the HES. Note that this incompleteness is intended, 
because the main aim of the search for metal-poor stars with 
the HES is to identify stars with [Fe/H] < -3.0. Therefore, the 
selection of candidate metal-poor stars was designed such that 
as many stars at [Fe/H] > -3.0 as possible are rejected, while 
main taining a high degre e of completeness at [Fe/H] < -3.0 
(see lChristlieb et al.ll2008l for details). 

For a star of a given [Fe/H], the corrections are also a func- 
tion of B - V color, being higher (more likely for a star to be 
included in the HES) for redder stars. The variation over the 
B - V color range of the HES sample can, in extreme cases at 
the higher metallicities, correspond to a variation of a factor of 
8 in selection efficiency (see, e.g., the line for [Fe/H] = -1.95 in 
Tab.©. 

For our comparison with the MDF of the Galactic GCs we 
adopt the [Fe/H] values from the current version of the on- 
line database of lHarrisI (1 19961) . The values for M15 and for 
NGC 7099 were updated with small corrections based on de- 
tailed abundance analyses carried out by J. Cohen and collab- 
orators (Cohen & Huang, in preparation; Cohen, Melendez & 
Huang, in preparation). The HES is (intentionally) incomplete 
for [Fe/H] > -2.0, so we only consider the set of GCs with 
[Fe/H] < -1.95, which contains only 16 clusters. We note that 
many analyses have shown that the Galactic GCs exhibit the 
same behaviour of abundance ratios (such as the increase of 
rC a/Fe] with decreas ing [Fe/H]) as the halo stars (e.g., Fig. 23 
of ICohen et al.|[2004h as do the halo stars. Thus, the conversion 
between a Ca line index and [Fe/H] adopted by the HES should 
be appropriate for Galactic GCs stars as well. 

Fig. [17] shows the cumulative MDF from the HES sample 
and for the Galactic GC system. The raw MDF and that cor- 
rected for the selection efficiency of the HES, given in Tab. [4] are 
shown. Note that the selection efficiency takes into account the 
uncertainties for [Fe/H] which result from the uncertainties of 
the measurement of KP and (B - V) from HES spectra, which re- 
sult in o"[ Fe / H ] between 0.2 dex and 1 .0 dex, depending on (B - V) 



Table 5. Cumulative halo MDF for [Fe/H] < -2.0 as observed 
in the HES (column "Raw"), and corrected for the selection effi- 
ciency of the survey (columns 3-5). For details see text. 
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colour (Christ lieb et al.l 120081) . Simulations suggest that for a 
sample with more accurate [Fe/H] determinations, such as the 
Galactic GCs, the maximum difference in the cumulative distri- 
bution up to [Fe/H] ~ -2.0 compared to that given in Tab. [4] 
is small and does not exceed the difference between the various 
corrected cumulative MDFs given in columns 3 to 5 of that table. 

The solid, middle line in Fig. [17] corresponds to the case 
where corrections according to the dereddened B - V color of 
each individual star of the HES sample have been applied. Since 
these corrections are themselves uncertain, two other variants are 
shown in this figure, and listed in Tab. [3] to indicate the poten- 
tial impact of the choice of B - V color on the corrections. The 
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Fig. 17. Cumulative MDF for [Fe/H] < -2.0 as observed in 
the HES (dashed line), and with corrections for the HES selec- 
tion efficiency applied (solid lines). Three different ways of ap- 
plying the corrections are shown to illustrate their uncertainty: 
Multiplication of the observed metal-poor star counts with the 
corrections for {B-V)o = 1.0 (upper solid line), the {B-V)o color 
appropriate for each individual star in the HES sample (thick, 
middle line), and the corrections for (B - V)o = 0.5 (lower solid 
line). The cumulative MDF of the GCs [Fe/H] < -1.95 is shown 
by filled circles. It agrees well with that of the halo field stars if 
the selection efficiency corrections are applied. 



first adopts the corrections for the bluest B - V color of Tab. |4j 
which are always the smallest, while the the second uses that of 
the reddest B - V color of Tab. [4] which are always the largest. 

Fig. [T71 shows that once the selection efficiency corrections 
given Tab.|4]are applied, the halo field star MDF we deduce here 
is a good match to that of the Galactic GCs. Instead of expecting 
roughly 10% of the sample covering the range [Fe/H] < -1.95 
to have [Fe/H] < -3.0, we expect only ~ 2 % to be this metal 
deficient, when the selection efficiency for the HES is taken into 
account. At [Fe/H] < -2.5, the expected fraction decreases from 
50 % to 8 %. Thus, the absence of any GC more metal-poor than 
-2.5 dex among a sample of 16 clusters at [Fe/H] < -1 .95 is not 
suprising. 

A similar situation holds for the stellar population in the 
dSph satellites of t he Galaxy. It has be en widely claimed (see, 
e.g. the review by iGeisler et al] 12007) that these dSph stellar 
populations show a significant la ck of stars with Fe- metallicity 
at [Fe/H] < -3.0. For example. iHelrni et al] d2006l) make this 
claim for the four systems for which they assembled the neces- 
sary data; i.e., Carina, Fornax, Sculptor and Sextans. 

Abundances are now available for large samples of stars in 
the nearest dSph galaxies. We concentrate here on those where 



there is little or no evidence for recent star formation and for 
which suitable samples are available. There are two additional 
issues that arise in a comparison of the stellar population of the 
dSph galaxies with the Galactic halo MDF. The first is that these 
metallicities are derived from line indices which measure the 
strength of the Ca infrared triplet (CaT) in moderate-resolution 
spectra. The conversion from a Ca abundance to a Fe abundance 
is a crucial issue, since the dSph stellar population clearly shows 
a different tren d of [Ca/Fe] versus J F e/H] than does the G alactic 
halo (see, e.g.. IGeisler et al]|2005l or lMonaco et al.ll2007l) . with 
[Ca/Fe] being smaller at a given Fe-metallicity in dSph galaxies 
as compared to GCs and the halo field. The second is how the 
sample to be observed spectroscopically in the dSph is selected. 
If, e.g., an equal number of stars in each color bin is chosen to 
probe the full range of color across the upper RGB in a dSph, the 
sample may be biased in metallicity, because the position of the 
upper RGB in the color-magnitude diagram depends on [Fe/H]. 
Instead, a representative subset of stars reflecting the color dis- 
tribution of the stars on the RGB should be chosen. 

Bearing these caveats in mind, we have constructed the cu- 
mulative MDF for several dSph galaxies. Given the larger sam- 
ples of very low-metallicity stars in these galaxies as com- 
pared to the limited number of such low metallicity Galactic 
GCs, we have normalized the dSph MDFs to [Fe/H] = -2.3. 
The selection efficiency of the HES over this lower metallicity 
range will be higher (i.e. closer to 1.0) and not vary as much 
than is the case over the regime we needed to consider for the 
Galactic GCs. The MDFs for Draco and for Ursa Minor (top 
row in Fig. [T8l are based on the database of IWinnickl (2003). 
She measured CaT line strengths from spectra obtained with 
the multi-fiber instrument Hydra at the WIYN telescope. Her 
sample is selected from radial-velocity members with no metal- 
licity bias. Winnick calibrates a relation between both [Ca/H] 
and [Fe/H] and CaT from observations of GC giants, making 
no attempt to take into account the difference in the behavior 
of [Ca/Fe] with [Fe/H] in these two stellar populations. We use 
her [Ca/H](CaT) values in t he figure (solid points) . The de- 
tailed abundance analyses of ICohen & Hu ang (2009) suggest 
that [Fe/H] ~ [Ca/H](CaT, Winnick) - 0.2 for Draco giants; 
this yields the open circles in the Draco panel of Fig. [18] The 
(constant) offset arises primarily from the lower [Ca/Fe] seen 
among dSph giants as compared to GC giant s , and i ncludes any 
error in the calibration adopted by Winnick| d2003l) . The MDF 
for Draco, with 24 stars at [Fe/H] < -2.3 (nine of which have 
[Ca/H] < -2.3), and for UMi, agree reasonably well with the 
halo Fe-MDF, once the selection efficiency of the HES is taken 
into consideration. 

For th e Fornax dSph we use the VLT/FL AMES +GIRAFFE 
survey of iBattaglia et al] d2006l) (their Tab. 4). The DART team 
in this paper converted their measurements of the strength of 
the infrared Ca triplet into Fe-metallicities (filled circles) us- 
ing the relation established by Rutledg e" et al] d 1997b . which was 
calibrated using globular cluster giants. Battaglia advises (priv. 
comm.) that their sample should be unbiased with respect to 
m etallicity. Although the re are only 7 Fornax stars in the sample 
of IBattaglia et al] (120061) with [Fe/H] < -2.3, the left panel of 
Fig- HH shows that the Fornax dSph is clearly deficient in such 
low metallicity stars relative to the MDF of the halo field stars 
when the HES raw counts are used. Once the selection efficien- 
cies are folded in, the Fornax cumulative MDF at [Fe/H] < -2.3 
is much closer to that of the Galactic halo field stars as inferred 
from the HES. 

The lower right panel of Fig. [l8]shows the cumulative MDF 
from the combined DART sample for the Carina, Sextans, and 
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Fig. 18. Cumulative MDF for [Fe/H]< -2.3 as observed in the HES (dashed line), and with corrections for the HES selection 
efficiency applied (solid lines; see the caption of Fig. [17] and the text for a detailed explanation). The filled circl es indicate the 
cumulative MDFs of Ursa Minor (upper left) and Draco (upper right), using the [Ca/H] values from Winnick (2003); for Draco we 
also show an estimate of the MDF for [Fe/H], assuming that [F e/H] = [Ca/H] - 0.2 (open circles). In the lower left panel we show 
a comparison of the cumulative MDF for Fornax from Batta glia et al.l (12006b . In the lower right panel the combined cumulative 
MDFs of Carina, Sextans, and Sculptor based on [Fe/H] values determined in the DART project ( filled circles; thi s data was kindly 
provided by the DART team) is compared to that of the HES. The result adopting the calibration of IWinnickl (120031) instead is shown 
as the open circles. 
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Sculptor dSph galaxies, with [Fe/H] values kindly supplied by 
the DART project (filled circles). This yields a total sample of 76 
stars with [Fe/H] < -2.3. Again, once the selection efficiencies 
are folded in, the cumulative MDF for these three dSph galaxies 
at [Fe/H] < -2.3 is much closer to that of the Galactic halo field 
stars as inferred from the HES than when the raw HES counts 
are used, but there still appears to be a deficit of stars in the 
combined dSph MDF at the lowest metallicity compared to the 
HE S MDF. 

Battaglia et al ] (120081) discussed the accuracy of their conver- 
sion between Ca triplet line strength and [Fe/H], given the differ- 
ence in the behavior of [Ca/Fe] with Fe-metallicity between GCs 
and dSph populations. Using a comparison of high-dispersion 
abundance analyses with their results from CaT measurements 
for a limited sample of dSph giants, they conclude that their Fe- 
metall icities are robust t o with in +0.2 dex. However, as pointed 
out bv lCohen & Huand (120091) . there are substantial differences 
betwee n the cali b ration adopted by the VLT DART project and 
that of IWinnickl (120031) . which suggest that the DART project 
metallicities are too high for [Fe/H] < -2.0. Hence, we con- 
verted the DART [Fe/H] values to those that w ould have bee n 
inferred using the CaT calibration to rFe/H] o f IWinnickl (|2003), 
combining Eq 13 of iBattaglia et all d2008l) with Eq. 3.5 of 
IWinnickl d2003l) . The results are indicat ed by the open c ircles. 
The application of the CaT calibration of Winnick (2003) to the 
DART data produces a better agreement with the HES Galactic 
halo MDF. Whil e the CaT techniqu e appears to be valid even at 
[Fe/H] < -2.5 (IStarkenburdl2009l) . the metallicity calibration 
needs to be improved in this [Fe/H] range. It is clear from the 
above that this issue is crucial in constructing a MDF. Efforts to 
validate and impro ve the calibration are currently underway by 
Starkenburg (2009) and others. 

We thus find that the MDF of the Galactic halo field stars, 
as derived from the HES, agree reasonably well with that of the 
Galactic globular cluster system and of the stellar population of 
the nearest dSph satellites of the Galaxy, when the calibration 
for converting infrared Ca triplet line strengths into [Fe/H] of 
IWinn ick (2003) is adopted. This holds over the range [Fe/H] < 
-2.3, after the selection efficiency corrections to the apparent 
MDF from the HES have been applied. However, ^- 2 -tests reveal 
that the differences between the halo MDF and the MDFs of the 
GC system and the dwarf satellites are still highly significant. 
If the original DART calibrations and [Fe/H] are valid, adding 
the HES selection efficiency corrections considerably improves 
the agreement in deduced MDF of the dSph galaxies with the 
Galactic halo field stars, but still leaves a problem at the lowest 
metallicities. 

Recently iKirbv et alj d2008l 120091) developed a spectral- 
synthesis technique that does not use the CaT at all. They found 
15 stars with [Fe/H] < -3.0 in seven of the ultra-faint dSph 
galaxies recently discovered by the SDSS. Since all these very 
low luminosity galaxies have mean [F e/H] values of - 1 .9 dex 
or lower, this is perhaps not surprising. ICohen & Huand f2009) 
have obtained high resolution spectra of a sample of stars in 
the Draco dSph, one of the more luminous of the dSph satel- 
lites of the Galaxy, and found one star with [Fe/H] < -3.0 in 
that dSph, in additi on to a Draco giant a t [Fe/H] = -2.97 dis- 
covered earlier by Shet rone et al.l (1 19981) . In Sculptor, one star 
with [Fe/H] < -3.5 has recently been identified based on high- 
resolution spectroscopy (Frebel 2009, priv. comm.). Finally, in 
a sample of 16 rad ial velocity members of the Bootes I dSph, 
iNorris et ail ([2008) have reported a giant with [Fe/H] = -3.4, 
based of measurements of the Ca II K line. Follow-up, high- 
resolution, high signal-to-noise observations with VLT/UVES 



confirm the result (Norris 2009, priv.comm.). Thus, extremely 
metal-poor stars are present, albeit in small numbers, in both the 
ultra-faint and classical dSph satellites of the Galaxy. 

8. Discussion and conclusions 

In Sect. [6] we have shown that a reasonable agreement with the 
overall shape of the HES MDF can be obtained for [Fe/H] > 
-3.6 by most models of Galactic chemical evolution, but only 
the model of Salvadori et al. with Z cr = 10~ 3 ' 4 Z Q reproduces 
the the sharp drop at [Fe/H] ~ -3.6 seen in the HES MDF. The 
lack of stars at [Fe/H] < -3.6 is highly significant: The mod- 
els typically predict that about ten such stars should be present 
in the HES sample, while only two are found. The significance 
of this discrepancy is reflected in the low probabilities for the 
MDFs predicted by the models and the HES MDF having the 
same parent distribution, as determined by KS-tests. It remains 
to be investigated whether the drop can be reproduced by modi- 
fying some of the assumptions of the models, or by adding fur- 
ther ingredients. 

The HES sample discussed in this paper contains no objects 
with [Fe/H] < -4.2, but considering the abundance analyses of 
three additional stars in this metallicity range published in the 
recent literature, it is obvious that it exists. However, a thorough 
and quantitative comparison with theoretical MDFs has to await 
larger statistically complete and unbiased samples which include 
more stars with [Fe/H] < -4.0. Such samples will become avail- 
able through new, deeper surveys for metal-poor stars that will 
comme nce in the near fu ture; in particular, the Southern Sky 
Survey (iKeller et al.l2007l) and a survey to be conducted with the 
Chinese 4 m Large sk y Area Multi-Ob ject fiber Spectroscopic 
Telescope (LAMOST; IZhao et al.ll2006l) . 

In the ACDM picture, the Galactic halo was largely built 
out of disrupted satellite galaxies. If stars had already formed 
within them at the time of accretion, then the MDF of the 
Galactic halo and of the existing dSph galaxies should agree 
at the metal-poor end with regard to the presence of a weak 
tail of stars with [Fe/H] < -3.0. It is thus encouraging for 
the ACDM scenario that our analysis shows better agreement 
between th e halo MDF and that of the dSph galaxies than 
claimed by iHelmi et al.l d2006l) . However, even if this were not 
the case, it would not necessarily be a strong contradiction to 
the ACDM scenar i o. Acc ordi ng to the semi-analyt i cal m odels 
of ISalvadori etail d2008l) and ISalvadori & Ferraral d2009l) . the 
MDFs of dSph galaxies can differ quite significantly from each 
other, depending on their individual enrichment histories. Hence 
their MDFs can also be different from that of the Galactic halo. 
An important question remaining to be answered is how the 
elemental-abundance ratios of the dSph stars at [Fe/H] < -3.0 
compare with those of the Galactic halo stars. 

Since the HES and the HK survey are in-situ surveys that 
predominantly sample the inner-halo population of the Galaxy 
(with R < 15 kpc), it is mandatory to consider the possibility that 
the (for now, poorly studied) outer-halo population of the Galaxy 
may indeed contain significant numbers of stars with [Fe/H] < 
-3.6, as might be indicated by the shift of the peak metallicity of 
the other-halo stars studied bv lCarollo et al.l (120071) to [Fe/H] = 
-2.2, a factor of four lower than the peak metallicity of inner- 
halo stars. This possibility is being actively pursued by high- 
resolution spectroscopic follow-up of stars that are likely to be 
members of the outer-halo population, based on their kinematics, 
by a number of groups. 
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